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Abstract 

a: 

CIh' We study production of photons via the non-equilibrium relaxation of a Dis- 



oriented Chiral Condensate with the chiral order parameter having a large 
initial amplitude along the 7r° direction. Assuming the validity of the low 



X 

energy coupling of the neutral pion to photons via the Ua(1) anomalous ver- 
tex, we find that for large initial amplitudes along the 7r° direction, photon 
production is enhanced by parametric amplification. These processes are non- 
perturbative with a large contribution during the non-equilibrium stages of 
the evolution and result in a distinct distribution of the produced photons 
and a polarization asymmetry. For initial amplitudes of the 7r° component of 
the order parameter between 200-400 MeV, corresponding to energy densities 
between 1-12 GeV/fm 3 we find a peak in the photon distribution at energies 



between « 300 — 600 MeV. We also find polarization asymmetries typically be- 
tween 5 — 10%. We discuss the potential experimental impact of these results 

as well as the problems associated with its detection. 
ll.30.Qc, 11.30.Rd, 11.40.Ha 



2 



Typeset using REVTeX 



I. INTRODUCTION 



The possibility that a disoriented chiral condensate (DCC) might form during the early 
stages of a high energy, heavy-ion collision, has attracted considerable attention in recent 
years PJ. The essential idea is that, near the central rapidity region, large energy densities 
~ a few GeV/fm 3 are achieved, corresponding to temperatures above 200 MeV, at which 
the chiral symmetry is restored f2|. As this region cools down via hydrodynamical expan- 
sion, it is conceivable that the QCD order parameter might be trapped in a configuration 
where it points in a direction different from the a direction which corresponds to the true 
vacuum. Subsequent relaxation of such a DCC configuration to the true vacuum is expected 
to be accompanied by a strong coherent burst of pions, which would then provide a defini- 
tive experimental signature of the chiral phase transition. The formation and relaxation of 
these DCC configurations could potentially provide an explanation for the CENTAURO and 
JACEE events found in cosmic ray emulsion experiments 0. 

Recent theoretical investigations |3|-|6| of effective field theories of QCD, particularly the 
linear sigma model seem to indicate that there is the possibility that such coherent regions 
may be formed during the non-equilibrium stages of the chiral phase transition. Results of 
classical semiclassical M and quantum calculations, with [H and without || expansion 
seem to indicate that there is a range of initial conditions for which the chiral phase transition 
goes through a non-equilibrium stage during which coherent pion clouds like DCC's may be 
formed. 

The telltale experimental signal for DCC events would be a non-binomial distribution for 
the ratio of the number of neutral pions to the total number of pions R = n n o / (n^o + n n + n -) 
with a probability distribution P(R) = 1/(2 V^) @. Both MINIMAX at the Fermilab 
Tevatron and WA98 at CERN are currently looking for strong fluctuations in the ratio of 
neutral to charged pions, and STAR and PHENIX at RHIC are scheduled to search for such 
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events. 

There are potential ambiguities in obtaining clear signals for the probability distribution 
P{R). Since the emitted pions are strongly interacting objects they will undergo strong final 
state interactions, and it could be that a DCC signal may be distorted and become indis- 
tinguishable from the background. An important limitation for the determination of P(R) 
is the accounting of neutral pions which decay to two photons before reaching the detector. 
The neutral pion distribution is typically reconstructed from the analysis of the emitted 
photons, a very demanding analysis in typical heavy-ion collisions with large multiplicity 
events. 

It is thus important to identify alternative signatures of the chiral phase transition, and 
DCC formation and relaxation. Electromagnetic probes such as photons and dileptons are 
ideal candidates since their interactions are weak, and hence they can carry the information 
from the early stages of the evolution of the plasma with minimal distortion. It is very likely 
that to extract exotic signals for the chiral phase transition or DCC events, both hadronic 
and electromagnetic probes will be needed. 

Furthermore, with the possibility of non-equilibrium effects during the formation and ex- 
pansion of the fireball in heavy ion collisions, the study of unusual non-equilibrium processes 
may provide a greater understanding of the signatures of the Quark-Gluon and Chiral phase 
transitions in experiments at RHIC and LHC. 

In an earlier work pj, we have recently studied direct photon production enhancement 
from the non-equilibrium relaxation of DCC's. There, we studied the gauged linear sigma 
model and concentrated on the electromagnetic interactions of the hadronic current. These 
calculations involved the electromagnetic interactions of the charged pions alone. The conclu- 
sion of this study was that the total number of photons produced during the non-equilibrium 
stages of a "quenched" phase transition (which lasts only for a few fms./c until the order 
parameter reaches its equilibrium value), can become comparable to the number of photons 
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released via (perturbative) decay of all the neutral pions produced in the heavy-ion collision 
(time scale for decay ~ 10 7 fm/c). 

In this article we focus our attention on the couplings of the neutral pion only, to elec- 
tromagnetism, through the axial anomaly. We work within the framework of the linear 
sigma model and introduce the anomalous interaction of the 7r° by means of the effective 
Wess-Zumino-Witten vertex. As described above, a DCC corresponds to a rotation of the 
chiral order parameter in a coherent domain. Since electromagnetism remains as an un- 
broken symmetry, such a rotation must be confined to the a — 7r° plane. A rotation that 
takes the order parameter out of the neutral meson plane homogeneously, will result in a 
global breaking of electromagnetism, which is not acceptable. Localized regions in which the 
charged pion fields obtain (locally) an expectation value (charged droplets) that average to 
zero on large volumes cannot be ruled out. However, the study of their dynamics involves 
a non-perturbative treatment of large amplitude inhomogeneous configurations, a very dif- 
ficult task and beyond the scope of this article. An initial step towards understanding the 
non-equilibrium evolution of small amplitude inhomogeneous configurations was presented 
in g. 

Minakata and Muller suggested in a recent work |J that a rotation in the neutral meson 
plane could indeed be induced in heavy-ion collisions. These authors argue that the presence 
of strong electromagnetic fields in these collisions will give a "kick" to the chiral order 
parameter in the 7r° direction, and that this effect is enhanced by the chiral U(l) anomaly. 
In this scenario, during the initial stages of the collision the order parameter acquires a 
component in the direction of the neutral pion. This is consistent with the formation of a 
"chirally disoriented" domain. During the expansion stage this disoriented region cools down 
and the order parameter relaxes to the true vacuum configuration, with the field oriented in 
the "a" direction. 

The goal of this article is to show that if one begins with the above-mentioned "disori- 
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ented" initial conditions with a large amplitude component in the ir° direction, the DCC 
relaxation will be accompanied by a coherent burst of photons which is released via a novel 
non-perturbative mechanism of parametric amplification, along with with the enhanced pion 
emission resulting from the non-equilibrium spinodal instabilities. 

The essential idea behind this novel mechanism is that the component of the homogeneous 
order parameter along the 7r° direction (the Fourier mode of the 7r° field with zero spatial 
momentum) couples to the electromagnetic field through the chiral anomaly, and effectively 
acts much as a time dependent "mass" term or "squeezing" term for the photon field that 
leads to the resonances in the equations of motion. 

This non-equilibrium mechanism exhibits the features of unstable bands, and exponen- 
tial growth of the photon mode functions, that are characteristic of parametric resonance. 
Furthermore, since the amplification and consequently the enhanced photoproduction will 
be confined only to certain frequency bands, this phenomenon translates into very peculiar 
distributions for the produced photons. Furthermore, because of the pseudoscalar nature 
of the interaction, there is a polarization asymmetry in the produced photons, i.e. photons 
with different polarizations are produced in different amounts. Both features of this non- 
perturbative mechanism might be potential signatures of DCC's or similar non-equilibrium 
relaxation phenomena during the chiral phase transition. Parametric amplification of pions 
has been conjectured to play a role in the hadronization stages of the QGP ||10|| . 

In this article we do not incorporate hydrodynamic expansion and treat the dynamics via 
the "quench" approximation HQ from an initial state in local thermodynamic equilibrium 
above T c . We also study the evolution of the system from a non-equilibrium initial state at 
Tj = to explicitly show that the qualitative features are robust and not tied to the "quench" 
scenario. Numerical results of the semiclassical evolution and quantum evolution || with 
expansion reveal that the qualitative features, such as pionic amplification and the estimate 
of the dynamical time scales are well described by the "quench" scenario. 
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Our aim in this article is to understand the qualitative aspects of this new mechanism of 
photoproduction in the simplest of settings, in light of the possibility that it could be a sig- 
nature of DCC's in heavy ion collisions. A deeper study of these and other photoproduction 
mechanisms including boost invariant expansion will be presented elsewhere [ p~5|j . 

This article is organized as follows: in Section II, we introduce the model and review the 
approximations used and their validity Section III is devoted to a simple derivation of the 
equations obeyed by the photon mode functions which in turn indicate the nature of the non- 
perturbative processes of anomalous photoproduction. In Section IV we review the numerical 
study of the equations of motion, and obtain the photon spectrum for a number of possible 
initial conditions, and in Section V we discuss some of the experimental implications of these 
results and potential difficulties for its detection. In Section VI we present our conclusions 
and future avenues of study. An appendix is devoted to a detailed derivation of the non- 
equilibrium Green's functions which are used for computing the expectation values of the 
relevant observables. 



II. THE MODEL 



In our subsequent investigations in this article we will work with the linear a model [|TT 
and furthermore, will take into consideration the electromagnetic interactions of the neutral 
pion only. The interactions of the charged pions give rise to out of equilibrium processes that 
produce a single (direct) photon at lowest order in a through spinodal growth of charged 
pion fluctuations; these processes have recently been studied in j7|. On the other hand, the 
anomalous interactions of the tt° will result in the distinctly different process of diphoton 
production, with the 2 emitted photons being polarized perpendicular to each other and will 
lead to a polarization asymmetry due to the pseudoscalar nature of the interaction. The 
Lagrangian density of the model under consideration is given by 
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(2.1) 



where $ is an 0(N + 1) vector, $ = (er, 7r°,7r) and 7? represents the iV — 1 pions, 7? 



(7T 1 , 7T 2 , 7T 3 , ....7T N ~ 1 ) and m 2 (t) introduces quench situation pi "by hand". As emphasized in 
|J the sudden quench approximation can be generalized to allow a "cooling" relaxation time 
(of the order of one fm/c) without major qualitative changes in the result and consistent 
with the numerical results including expansion 
We then take 



We will also consider below the evolution of the system from a non-equilibrium initial 
condition at T« = for comparison. 

The last term in eq. (|2.1|) results from the axial anomaly. In a constituent quark model 
it is known to arise from the quark triangle diagram, which for 3 flavors of quarks leads to 
the amplitude for 7r° decay into two photons described by the effective interaction vertex in 
the above Lagrangian ||12|| . 

At this point we want to remark that while the nature of the effective Lagrangian for the 
neutral pion decay into photons is clearly understood in perturbation theory, it is by no means 
clear that this effective interaction also describes strongly out of equilibrium situations. It 
has been recently pointed out by Pisarski, Tytgat and Baier et. al. [[L3[] that in equilibrium 
and above the critical temperature there is a drastic modification of the anomalous vertices. 
It is conceivable that the anomalous decay vertices are modified in strongly out of equilibrium 
situations as well. We are currently studying this possibility in detail within the framework 
of non-equilibrium field theory in a constituent quark model, but in this article we accept 
the anomalous vertex determined by the perturbative triangle diagram in equilibrium and 
study its consequences. 



m 2 (t) 




Q(-t) - 1 Tj > T c . 



(2.2) 
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The linear sigma model is at best a phenomenological model, and its parameters are fixed 



by the low energy phenomenology of pion physics to be |TT| : 

M a « 600 MeV ; / w » 93 MeV ; A « 4.5 

/i » (120 MeV) 3 ; T c « 200 MeV. (2.3) 

This value of the critical temperature T c is a consequence of the parameters of the model 
and is somewhat larger than the lattice estimates T c « 150 MeV. This is a cutoff theory with 
an ultraviolet cutoff of the order of A « 1 GeV. Although we will be dealing with the 0(4) 
model with 3 pions, we will use the large N limit to provide a consistent, non-perturbative 
framework to study the dynamics. For a clear description of this approximation and its 
validity, and its applications to a wide variety of non-equilibrium situations we refer the 



reader to [ I^ , [TB|JT7|| . Eventually one should set A^ = 3 which is the number of physical 
pions in the real world and though questions may be raised as to the validity of the large 
A" expansion in this case, it provides a systematic approximation method to study the non- 
perturbative physics of the sigma model. More importantly perhaps, it is an approximation 
scheme that preserves energy conservation (in the absence of expansion), implements PCAC 
and satisfies the Ward identities of chiral symmetry. As a result, the pions in the asymptotic 
equilibrium state are exactly massless in the absence of the explicit symmetry breaking term. 

Non-equilibrium quantum field theory requires a path integral representation along a 
complex countour in time |IS| : 19]| , with the Lagrangian density along this contour given by 



C neq = £[$+] -£[$-] (2.4) 
with the fields $ ± defined along the forward (+) and backward (— ) time branches. For 



further details see references [14,l7|. The non-equilibrium equations of motion are obtained 



via the tadpole method |I7J. As mentioned earlier, in the situations of interest to us in this 
article, both the a and 7r° fields may acquire expectation values. Accordingly, we first shift 
a and tt° by their expectation values in the non-equilibrium state 
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a(x, t) = <P(t) + X (x , t) ; 0(t) = (a(x, t)) (2.5) 

7r°fo*) = C(f)+#M) ; at) = (n°(x,t)) 

with the tadpole conditions 

( X (x,t)) = ; (^(x,t)) = ; (vf(f,t)) = (2.6) 

being imposed to all orders in the corresponding expansion. To leading order in 1/N, the 
large-iV approximation is equivalent to the following Hartree factorization of the quantum 
fields: 

X 4 , (i/) ) — > 6(x 2 )x 2 + constant, (6(ip 2 )ip 2 + constant) 

X 3 , (^ 3 ) - 3( X 2 )X, (3<V», (2.7) 
4 

(tt • tt) 2 -> (2 + — )(7r 2 )7f 2 + constant 

TfV, (*V) - 7T 2 ( X 2 ) + (7r 2 ) X 2 , (t? 2 ^ 2 ) + <vr> 2 ) 

7f 2 X , (ttV) - (7r 2 ) X , ((vr». (2.8) 

In addition to the above factorizations, we will implement a further Hartree factorization of 
the Wess-Zumino term, 

i)FF ^ip(FF). (2.9) 

There is no a priori justification for the factorization of the anomalous term, but this allows 
us to treat the full back-reaction problem self-consistently maintaining energy conservation 
(in the absence of expansion). 

To leading order in 1/N, the effective non-equilibrium Lagrangian density is then given 

by 

£[0 + X \ 7T+] - £[0 + x", tt-] = {\(d X + ? + \(d^ + f + \(d^f 

-Vi(t) x + - v 2 (t)^ + + K{tm)x4> - l -M x +2 {t) X +2 - ^av 2 w +2 - l -Mt +2 {t)t+ 2 } 

-\f; v F^ + ^C( t )F+F + '" + ^i> + (F+F + n -{+—>-}. (2.10) 
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where 

Vi(t) = 4>(t) + 0(t)[m 2 (t) + 4A0 2 (t) + 4AC 2 (t) + £(*)] - /i (2.11) 



V 2 (t) = C(t) + CW[m 2 (t) + 4A0 2 (t) + 4AC 2 (t) + £(£)] 



(2.12) 



A<*(*) 



m 2 (t) + 4A0 2 (t) + 4AC 2 (t) + £(t) 



(2.13) 



m 2 (t) + 12A0 2 (t) + 4AC 2 (t) + E(*). 



(2.14) 



Ml(t) 



m\t) + 12AC 2 (t) + 4A0 2 (t) + E(t). 



(2.15) 



E(f)=4A (7f 2 )(t)-(7r 2 )(0) 



(2.16) 



The subtraction of (7r 2 )(t) at i = leads to a renormalization of the mass term given in ( |2.2| ) 



. The presence of the mixing term does not contribute to the dynamics of the theory 
in the large N limit because such a term is formally of 0(1/ N) and is therefore subleading. 
In the large N limit, the Lagrangian becomes quadratic at the expense of a self-consistency 
condition for 

The pion fields obey a linear Heisenberg evolution equation (in terms of the self-consistent 
field), and the Heisenberg field operators can be expanded as 



where the destruction and creation operators of Fock states associated with the 

pion field and Q is the quantization volume. The frequencies W^ i determine the initial state 
and will be specified below. 

It is found mi that the pion mode functions £/&(£) and the order parameter(s) <j>(t), £(i) 
obey the following equations of motion 




[a k U k {t)e ik - s +a\Ul{t)e- ik - s \ 



(2.17) 
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4>(t) + [m 2 (t) + 4A0 2 (t) + 4AC 2 (t) + X(t)]<j>(t) -h = 



(2.18) 



C(t) + [m 2 (t) + 4A0 2 (t) + 4AC 2 (t) + E(f)]C(*) - ^(FF) = O (2.19) 



^ + k 2 + m 2 (t) + 4A0 2 (t) + 4AC 2 (t) + £(f) 



= (2.20) 



At this stage we must specify the initial conditions on the quantum pion fields. We will 
study the evolution with different sets of initial conditions to determine the sensitivity of the 
qualitative features of the dynamics to changes in the initial conditions: 

1. In the first case we will assume that at the initial time (the time of the "quench") the 
quantum fluctuations are in local thermodynamic equilibrium at the initial temperature 
Tj > T c and that the mode functions £4(0) correspond to the instantaneous positive 
frequency modes with mass 

M 2 = m 2 (t < 0) + 4A(0 2 (O) + C 2 (0)). (2.21) 

In addition, the order parameter is assumed to have a non-zero initial amplitude in the 
7T° direction i.e. C(0) ^ and 0(0) = 0. 

2. In the second scenario we assume that the order parameter is initially displaced from 
the equilibrium position exactly as in (1), but the modes are chosen to be at zero 
temperature, Tj = 0, so that the m 2 (t < 0) = — ^f-, and the mass for the modes is 

M 2 

M 2 = — + 4A(0 2 (O) + C 2 (0)). (2.22) 

However, in both the above cases, the modes must satisfy, 

U k (0) = 1; U k {0) = -iW k;i ; W Ki = W + W. (2.23) 
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These initial conditions are to be treated only as assumptions that allow concrete calcula- 
tions and estimates and are by no means a definitive description of the physics. The correct 
set of initial conditions should in principle be obtained by evolving parton cascade mod- 
els well into the regime in which hadronic physics can be described by low energy effective 
lagrangians. This is certainly an ambitious task and clearly beyond our current understand- 
ing. Thus in view of the lack of clear and precise knowledge of the initial conditions, we are 
restricted to investigating a physically plausible initialization of the dynamics. The above 
choice is the simplest that allows a concrete calculation and a quantitative description of the 
dynamics. 

With these initial conditions we are led to the self-consistent equation [|]|| 

M ■ (2M) 

For the initial conditions (2), we must set Tj = 0. In addition, the initial conditions on 
the order parameters and their derivatives, namely 0(0), 0(0), £(0) and £(0) must also be 
specified. At the moment they are completely arbitrary, constrained only by the requirement 
that typical energy densities in a DCC region are ~ a few GeV/fm 3 (which is the typical 
energy deposited in the central rapidity region in a heavy-ion collision). 

Clearly, a self consistent solution to the above equations will only be possible after un- 
derstanding the time evolution of the (FF) term in ( |2.19| ). This in turn, will require an 
investigation of the (FF interaction in ( |2.10| ). As shown below, it will give rise to out of 
equilibrium photon production via a process of parametric resonance. 



III. PHOTON PRODUCTION 

Let us now turn to the electromagnetic sector of the theory. Since the processes of interest 
to us do not involve photons in intermediate states, we can avoid potential gauge ambiguities 
by working in the Coulomb gauge and restricting our attention to the transverse, physical 
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degrees of freedom. In this gauge, the quadratic part of the Lagrangian for the physical 
photons, is given by 

^ = \VM* + ^fFF. (3.1) 

Since 

FF = 8(9 A)(6» j A fc )e°^ fc + total divergence (3.2) 

C q em may be re-written as, 

£L= \\dM 2 + ^^(doA^A^ (3.3) 

= — A T ■ UA T - —^lALd j M.e ijk + total divergence 
2 4vr 2 2/ w T 3 T 

leading to the Heisenberg operator-equation of motion: 

UA T + x At — 0. (3.4) 

Recall that the interaction term ipFF has been Hartee-factorized to give tp{FF) and does 
not contribute to the operator equation of motion above. 

The physical content of the equations of motion becomes manifest once we expand the 
transverse photon field in terms of the circularly polarized states. For a mode with wave 
vector k, consider the triad of vectors k, ei(fc), e 2 ik) with the properties 

ei(-fc) = -ei(fc) ; e 2 (-fc) = ? 2 (fc) (3.5) 
k x ei(fc) = e 2 (fc) ; k x e 2 (/c) = —ei{k). 

We can now construct the circular polarization vectors 

In terms of these circular polarization vectors we can perform the usual mode expansion 
for the fields in terms of annihilation and creation operators: 
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A T ^t)=^=^A T (k,t) (3.7) 

= ^ E Ty= { [M^is(*)e + $ + M^OM*)] e** + h.c.} 
We find that the mode functions obey the following equations of motion 

^ + ev u ( t )- k ^fv^0 (3.8) 
^ + ^( i) + ^f^0. (3.9) 

To these equations of motion we must append initial conditions. Now, we want to 
describe the process of emission of physical, circularly polarized massless photons. Therefore 
the initial conditions for these mode functions must be: 

^i,2fc(0) = 1 , \>i, 2fe (0) = -ik. (3.10) 

Furthermore, we will assume that at the initial time, the photons are in local thermody- 
namic equilibrium at the initial temperature. We will also study the case in which there are 
no initial photons by taking the initial temperature of the photon distribution function to 
zero. 

There are several noteworthy features of these equations of motion that deserve comment: 

1. The zero mode of the neutral pion field acts as time dependent mass term (with a 
momentum dependence), which will result in a "squeezing" of the quantum fields. 

2. The equations decouple when written in terms of the circular polarization mode func- 
tions. 

3. The effective "mass terms" for the two polarizations have opposite signs. These last 
two features are a direct consequence of the pseudoscalar nature of the coupling of the 
neutral pion to the vector field. 
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As the 7T° zero mode executes quasiperiodic oscillations, the solutions to the mode equa- 
tions ( |3.8| , |3.9| ) will display the phenomenon of parametric resonance, for some values of k. In 
other words, for certain (forbidden) bands of values in /c-space, the modes will grow exponen- 
tially with Floquet-like solutions. This growth of modes will be reflected in the exponential 
growth of particle number (in this case, the number of asymptotic photons). Furthermore, 
since the coupling affects right and left handed polarizations differently, there will be a po- 
larization asymmetry in the photon emission. A quantitative analysis of this phenomenon 
will require a solution of the coupled differential equations (|2.18|) , (|2.19| ), fl2.20| ), ( |2.24|) and 



( ft .8|, p.9|) which will be done numerically below. 

The number operator for asymptotic photons is given by 



A T (k, t) ■ A T (-k, t) + k 2 A T (k, t) ■ A T (-k, t) 



(3.11) 



In order to obtain the expectation value of this Heisenberg operator in the initial density ma- 
trix, one needs the quantum two-point functions for the theory. These can be obtained from 
the expansion (|3.8|) . A detailed derivation of these correlators in terms of non-equilibrium 
Green's functions has been given in an appendix. From the expressions for the Wightman 
functions ( |A18| ) and ( [A 19] ) it is easy to see that the expectation value of the number operator 



in terms of the mode functions (|3.8j , |3.9| ), has the following form, 

k 



(W)) 



1 
1 



coth 
coth 



2Ti 
k 



271 

= N+(k,t) + N_(k,t) 



\V lk (t)\ 2 + k 2 \V lk (t)\< 



V 2k (t)\ 2 + k 2 \V 2k (t)f 



1 

2 
1 
2 



(3.12) 



Only the symmetric part of the Green's functions contribute to the expectation value. From 
this expression for the number of asymptotic, transverse, massless photons produced, and the 



evolution equations (|3.8|J3.9| ) it is clear that the number of left and right circularly polarized 
photons produced by the evolution of the zero mode of the neutral pion field will in general 
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be different, since the evolution of the modes are different. The polarization asymmetry in 
the photons produced may be an experimentally relevant signal. We define this asymmetry 



as 



E(k,t) 



(N + {k,t) - N.(k,t)) 



(3.13) 



(N+(k,t) + N.(k,t)Y 

The dynamics of the process of photon production via the non-equilibrium evolution of 
the expectation value of the neutral pion field is then obtained by solving the self-consistent 
set of equations of motion (|218U2.19|JE2G|jE2l |3T8|j3Tg|) with the initial conditions (ET23, KW). 



IV. NUMERICAL ANALYSIS 

It is convenient to work with dimensionless variables for purposes of numerical analysis. 
The relevant dimensionful scale in this problem is Aqcd ~ 200 MeV = Mf, and we use it 
to scale out the dimensions of all the relevant variables in the equations of motion: 



/ x C(t) (f) k „ M 2 , , 



U q {t) = ^, W ( ) = — - , -^L(0) = -i^j- (4.2) 

xiMt) = , xiM Q ) = y=> = -l v^ ( 4 - 3 ) 

M| Mp Mp 

Thus the Wronskian condition for the mode functions defined in terms of dimensionless 
variables is 

-fr-Xito - -^rXiA ~ "2«. (4-5) 
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The expectation value of the anomaly term which appears in the equations of motion can 
now be obtained as follows 



(FF) = 8(d A^ r d j A k T )e ijh 



(4.6) 



8- 



d 



d 3 k r d 3 k 



^ +P) "(^(fc,t)^(fc',t'))(-^)e^|t=t'. 



dtJ (2n) 3 J (2tt) 3 

The two point correlator in the above expression can be obtained from the coincidence limit 
of the Green's functions of the full quantum theory, given in eqs. ( |A18j ) and ( A19 ). In this 
case, only the antisymmetric part of the Green's function contributes to the expectation 
value. The expectation of FF in the initial density matrix, is then found to be: 



(FF) 



Ml 
ir 2 



q 3 dq coth 



2r 



d 

dr 



(\X2q\ 



\Xlq\ 



(4.7) 



where q is the dimensionless momentum and use was made of the Wronskian condition 

The full dynamics of the system can now be followed by numerically solving the following 
set of equations: 



d 2 r](r) m 2 (r) 
~dP r + ^Ff 



4X V (T)[r ] 2 (t)+e(r) + J:(r)]-H = 0, 



(4i 



d 2 gr) m\t) 
dr 2 M 2 F 



£(r)+4A£(r) ^ (r)+ £ 2(r) + E(r) _ r(r) = 



d 2 U, 

~dV 2 



" + ( q 2 + ^ + 4A^ 2 (t) + 4A£ 2 (i) + 4X2(t)) U q (t) = 



M 



d 2 X i q , ( 2 aM F d£\ 

1 ~ -^^<1— I Xlq = 



dr 2 



d 2 X2g 
dr 2 



' 2 aM F di\ 
, n U dr J 
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with the initial conditions 

U q (0) 



:, W,(0) 



(4.9) 



Xl,2g(0) = — , Xl,2g(0) 

and the self-consistent fluctuations given by 



4vr 2 



/ 5 g 2 dg(|W g (r)| 2 -|W 9 (0)| 2 )coth 
Jo 



L0 n 



2r 



(4.10) 



r(r) 



q 3 dq coth 



2r 



a /. . 



8tt 3 ^ y 

The number of photons per unit volume per unit phase space is then 



(27T) 



,d({%))-{JV,(0)))/n_n( g) r) 



1 



4g 



coth 



Q_ 

2r 



d 3 q 
dx 



\Q\ 



1 7 1 2 _|_ | ^X2q 1 2 



+ q 2 (\Xi q ? + \X2 q \ 2 )-M 



(4.11) 



Here, we have taken the continuum limit, and subtracted out the number of photons present 
at t — which forms a thermal background at the initial temperature Tj. We have also 
introduced the quantity n(q,r) = (2tt) 3 \q\ ^Mllzi^YiMli/^ ^ Our numerical results for the 
photon spectrum and the invariant rate will be expressed in terms of this quantity. The 
invariant rate for photon production can be obtained by looking at the time derivative of 
( |4.11| ) and using the mode equations, it is found to be: 



(2vr 



(4.12) 



A comparison of this expression for the rate and that for the expectation value of the anomaly 
term in eq.( |4.7| ) clearly indicates the the important role played by the latter quantity in 
ensuring energy conservation. The photons produced lead to a non-zero expectation value 
for FF which then feeds back into the equation of motion for the zero mode. 
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The above expression for the rate will also turn out to be useful for obtaining quantitative 
estimates of the location of the unstable bands. In equilibrium, |xi<j| 2 = \X2q\ 2 = 1/?, and 
thus the rate vanishes. However, out of equilibrium and particularly for the processes that 
we have considered here, both \xiq\ and \x2q\ will acquire small time dependent modulations 
around a mean value, with a frequency that is roughly ~ q, since that is the relevant scale in 
the problem. Thus the oscillation frequency of J^dxi^l 2 — |X2q| 2 ) is ~ 2g. Now the maximum 
buildup of particles will occur when this oscillation is in phase with the derivative of the 
zero mode £(t), resulting in a positive rate for all times, causing a resonant production of 
photons. Therefore the unstable bands will be determined by the condition 

q res M F = k res « (4.13) 

with u n o being the oscillation frequency of the expectation value of the 7r° component of 
the order parameter. For small oscillations in the 7r° direction, u) n o = M n = 140 MeV and 
so the resonant band should be at k = 70 MeV, which is exactly the energy of two back 
to back photons produced by a tt° at rest. The validity of these qualitative arguments and 
particularly that of eq. (|4.13| ) will be borne out by the results of the numerical analysis 
which we present below. 

For all the cases that we considered, we chose to represent the "quench" from an initial 
temperature Tj = 1.1T C = 220 MeV to zero temperature, and we tracked the evolution of the 
system up to a time of about 10 fm/c which is the typical time scale after which the effects 
of hydrodynamical expansion become important. Since we are interested in studying the 
optimal situation for large anomalous photon production from the non-equilibrium evolution 
of the neutral pion, we have performed our numerical analysis in all cases with the initial 
value of the a component as well as its initial time derivative being zero. An initial condition 
with both the neutral pion and a components being different from zero will result in a smaller 
photon yield. This is because, for a given initial energy density, if the a component is large, 
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the 7T° component must be smaller and it is the 7r° zero mode that drives the parametric 
amplification. Furthermore since the sigma model is an effective theory below 1 GeV we 
look only at photons with momenta less than this value. Photons with energy larger than a 
GeV will be indinstinguishable from the hard photon background and much less significant 
as a signature of low pr physics associated with DCC's or the chiral phase transition. 
Results: 

The results of the numerical calculations are shown in figures (1) through (6). We first 
look at the result obtained for small oscillations in the 7r° direction with initial conditions 
((0) = £(0)M F = 0.1M F , f(0) = 0, 0(0) = 0(0) = 0. This is the perturbative regime, 
in which £(r) undergoes small amplitude (almost undamped) oscillations, with emission 
of back to back photons at 70 MeV. Fig.(l) shows the spectrum of photons emitted per 
unit volume, n(q,r) = (27r) 3 |g| d ( N ( T )-^(°))/ n ; (j n units of fm _1 ) at time t = 10 fm/c, for 
the small amplitude initial conditions described above. Here, we have subtracted out the 
thermal distribution of photons that are present initially at t — when the system is in 
equilibrium at a temperaure, Tj = 220 MeV. Clearly, there is a resonant peak which upon 
closer inspection is found to be exactly centered around 70 MeV. This corresponds to the 
usual perturbative decay of zero momentum pions and as expected, the peak value of the 
photon number, 2.75 x 10~ 5 , is 0(a 2 ). In this perturbative (small amplitude) regime the 
location of the unstable band for small oscillations can be easily obtained by solving the 
mode equations ( |3.8| , |3.9|) assuming a sinusoidal time dependence for the zero mode, which 



leads to a Mathieu equation for the non-zero modes for which the position of the unstable 



bands can be found in textbooks 20 



Figs.(2.a) and (2.b) show £(r) and n(q,r) = (27r) 3 |g| d(jV(r) ~^ (0))/ " , (in units of fni^ 1 ) at 
time t — 10 fm/c respectively, for the initial conditions £(0) = £(0)M^ = l.OMp, £(0) = 0, 
0(0) = 0(0) = 0. This corresponds to preparing a configuration with an energy density 
A(C(0) 2 — f 2 ) 2 ~ 0.5 GeV/fm 3 above the thermal energy density. £ now evolves with a larger 
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amplitude and also a higher frequency as shown in Figure (2. a). The small damping has its 
origins in mode mixing and excitations of the quantum fluctuations (pions and "sigmas") 
via spinodal and parametric amplification, and also to a smaller extent, the contribution 
from photon production. The spectrum of emitted photons is shown in Fig.(2.b). It features 
a prominent peak centered around 306 MeV with a much larger value, 0.005 fm _1 (vs. 
0.0000275 fm _1 in Fig.(l)) which is of 0(a) (rather than 0(a 2 ) as in Fig.(l)) showing 
explicitly the non-perturbative nature of the amplification phenomenon. The band is also 
quite broad with a width « 120 MeV. This is clearly a non-perturbative effect caused by 
parametric resonance. In fact our qualitative argument ( |4.13|) for obtaining the resonant 
frequency works rather well here. u n o can be read off from the oscillations of the zero mode 
in Figure (2.a) and is « \ x 200 MeV so that k res « « 310 MeV. To illustrate the 

point of parametric amplification, Fig.(2.c) shows the invariant rate h(q, r) = (2Tc) 3 \q\ dN j l PJ n , 
(in units of fm -2 ) for q at the center of the resonant band k = k res as a function of time 
(in units of fm/c). It is important to remark that unlike the small amplitude regime, in 
this large amplitude case, the rate is 0(a) rather than the perturbative result 0(a 2 ). This 
figure clearly displays a positive quasiperiodic solution with a growing envelope (the "true" 
invariant rate should be obtained by dividing by the total number of photons as a function 
of time, as the above rate has been defined as the time derivative of the number of photons, 
not of its logarithm). Our numerical studies show that, away from the center of the resonant 
band, the rate samples negative values also, thus leading to a lower photon yield for q values 
away from the center. 

With a larger initial amplitude for the 7r° zero mode, we see the same non-perturbative 
effects as before with a much larger enhancement. Figs. (3. a), (3.b) and (3.c) show the 
results for initial conditions C(0) = f(0)M F = 2.0M F , f (0) = 0, 0(0) = 0(0) = 0. The 
energy density corresponding to this initial state is ~ 12 GeV/fm 3 . It is not surprising that 
for such a large energy density, the peak in the photon spectrum (Fig.(3.b)) and the invariant 
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rate at the resonant frequency (Fig.(3.c)) are about an order of magnitude larger than the 
previous case. Furthermore the position of the peak is now shifted further to the right to 
625 MeV in qualitative agreement with the argument ( [4.131 ). 

As anticipated in the discussions of the previous section, these non-perturbative effects 
and the pseudoscalar nature of the interaction, will give rise to a polarization asymmetry in 
the photons produced, as well as a non-equilibrium expectation value for FF. 

In Figures (4. a) and (4.b) we show the polarization asymmetry at the resonant value 
of q; E(q res ,r) = (n + (q res ,r) - n_(q res , r)) / (n + (q res , r) + n_(q res , r)) as a function of time 
for k = 306 MeV, £(0) = 1.0 (corresponding to the initial conditions of Figs. (2)) and 
k = 625 MeV, £(0) = 2.0 (corresponding to Figs. (3)). Although this asymmetry is a (quasi) 
oscillatory function of time, the average over a time scale of about 10 fm/c, yields an estimate 
of |S(g res ,i)| = 0.05 — 0.10 for the values used in the figures. 

The behaviour of T(r) = aM F (FF) /8irf n (in units of M F ) as a function of time up to 
10 fm/c is shown in Figures (5. a) and (5.b). Though the behaviour is oscillatory, it is not 
symmetric about zero and is of order 10~ 4 — 10 -3 Mp. Thus (FF) itself is of order 0.1 — 1 
M F . 

We have also studied non-equilibrium initial conditions that do not correspond to lo- 
cal thermodynamic equilibrium by setting Tj to zero, in the initial conditions and relevant 
quantities. 

Fig. (6. a) shows n(q,r) = (27r) 3 |g| d ( JV (' r )~^( ))/ n ; (j n units of fin" 1 ) at time t = 10 fm/c 
respectively, for the initial conditions C(0) = £(0)M F = l.OMp, f (0) = 0, 0(0) = 0(0) = 0, 
Ti = 0. Comparing this figure with Fig.(2.b) we see that the main features are qualitatively 
the same. The amplitude is smaller in Fig. (6. a), reflecting the Bose-enhancement for 7^ 
in Fig.(2.b). There is also a slight shift in the center of the unstable band. In Fig. (6. a) the 
peak is at k = 284MeV rather than at k — 306MeV as in Fig.(2.b). Fig.(6.b) shows the 
invariant rate h(q, r) = (2TT) 3 \q\ dN ^J n , (in units of fm -2 ) for q at the center of the resonant 
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band k = k res as a function of time (in units of fm/c) for the same initial conditions as for 
Fig. (6. a). Fig.(6.c) shows the polarization asymmetry for the same initial conditions. 

Clearly the effect is fairly robust in that it mainly depends on the non-equilibrium aspects 
of the evolution with a rather small quantitative change due to Bose-enhancement if the 
initial state is in LTE at Tj 7^ 0. A similar conclusion is obtained after numerical analysis 
with the initial conditions of Figs.(3.a-c) and Tj = 0, with the same qualitative features 
modified slightly by Bose-enhancement (or the lack thereof). 

V. POSSIBLE EXPERIMENTAL SIGNATURES: 

In a heavy-ion collision one of the most important contributions to the production of 
low energy photons is from the anomalous single particle decay of the neutral pion. The 
neutral pions formed during the hadronization and chiral symmetry breaking stages will 
mainly decay at freeze out, since their lifetimes (at rest) are ~ 10 7 fm/c. 

The mechanisms studied above are effective only during the typical time scales for the 
non-equilibrium effects and for large initial amplitudes of the neutral pion component of 
the order parameter. Numerical estimates of the expansion time scales in the semiclassical 
Q and quantum large N approximations within the sigma model H reveal a longitudinal 
expansion time scale between 5 to 10 fm/c (spherical expansion scales are somewhat shorter). 
Although we have not incorporated expansion in our studies, it is a plausible assumption 
that even if the initial amplitude for the neutral pion component is large, giving rise to 
the non-perturbative parametric amplification phenomena described above, this amplitude 
will have diminished substantially by the end of the expansion stage. Thus most of the 
non-equilibrium, non-perturbative phenomena will take place during the initial stages for a 
time-scale of about 10 fm/c. Only the photons produced during this non-equilibrium stage 
will feature the abnormal distributions studied above, whereas the photons emitted via the 
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perturbative (small amplitude) decay of the neutral pion after freeze out will have the typical 
distribution with a peak at 70 MeV for decay at rest. 

Therefore an experimentally meaningful quantity is the ratio of the total number of pho- 
tons produced during the non-equilibrium stages with abnormal (parametrically amplified) 
distributions to half the initial number of 7r°'s. Notice that this assumes that all the 7r°'s 
produced will decay perturbatively to 2 photons. Thus this ratio gives a conservative esti- 
mate of the experimental relevance of the abnormal electromagnetic signal and provides an 
estimate for the experimental signature of the non-equilibrium aspects of anomalous photon 
production during the chiral phase transition. 

We found numerically that the pion number density produced via parametric ampli- 
fication and spinodal instabilities during the non-equilibrium stages, including the pions 
present in the thermal bath is about 0.6/fm 3 for ((0) = £(0)Mp = l.OMp; and ~ 5/fm 3 for 
((0) = £(0)M F = 2.0M F , where M F = 200 MeV. The corresponding numbers for the pho- 
tons produced by parametric amplification out of equilibrium are 0.0015/fm 3 and 0.005/fm 3 
respectively. 

Thus the ratio of the number of amplified photons to the number of decay photons is 
~ 0.1 — 0.2%, a rather small number. If the detector can measure final state polarization 
of photons, this signal must be correlated with the total (integrated over the lifetime of the 
"fireball") polarization asymmetry E tot = \n+ — + n_) ps 0.05 — 0.1 where the sum 

in the denominator refers to the parametrically amplified photons. 

Therefore the experimental signature of these non-equilibrium effects resulting in an 
enhancement of photons produced through the non-equilibrium decay of the neutral pion, 
will be a peak at momenta k > 300 MeV with an amplitude less than a percent of the 
amplitude of the usual n° peak and with a polarization asymmetry between 5 — 10%. (Recall, 
however that we have assumed that all the neutral pions produced will contribute to the 70 
MeV peak by decaying at rest, to two photons.) The photon energy corresponding to the 
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center of the peak increases with the initial energy density at the onset of the non-equilibrium 
stage. 

Clearly these are small effects leading to less than 1 percent signal over background 
relative to the usual n° peak. However, our estimate is only a lower bound and the effect 
could still be within the limits of resolution of the detectors and could potentially provide 
one of the electromagnetic signatures of DCC's. 

In the energy region k > 280 MeV (corresponding to an initial value of (ir°) 200 MeV) 
the most problematic interference for detection is the decay of the rj pseudoscalar meson into 
two photons, with a branching ratio of about 40%. This neutral mode of decay for the r\ 
meson has the same features as that of the neutral pion, insofar as angular distribution and 
polarization asymmetry are concerned. We can provide a rough (albeit naive) estimate of the 
number of diphotons produced by r\ decay by assuming that the 77's are in equilibrium at the 
initial temperature since they would not undergo parametric or spinodal amplification; this 
yields an 77 number density of about 0.088/fm 3 . This in turn leads to about 0.035/fm 3 di- 
photons at an energy w 275 MeV. Thus if the initial condition corresponds to a disorientation 
in the 7r° direction with an expectation value of the neutral pion field ps 200 MeV, this would 
translate to a peak in the photon spectrum in the vicinity of the peak associated with the r\ 
meson with about 5% ratio of the peak values. Clearly this is not a generic initial condition; 
for larger values of the 7r° component, the peaks in the distribution will move to larger values 
of the momenta and their amplitude will be enhanced. Although one could argue that for 
larger values of the momenta the background photon production from higher mass resonance 
decays will be important, these will be suppressed by thermal factors and smaller branching 
ratios. 
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VI. CONCLUSIONS AND FURTHER QUESTIONS 



In this article we have focussed on the study of photon production enhancement through 
the non-equilibrium stages of relaxation of a DCC via the C/a(1) anomaly. 

The premise of the study is that if during a heavy-ion collision a state of large energy 
density is formed in which the chiral order parameter is "disoriented" with a large amplitude 
in the 7T° direction, such a state will relax via parametric and spinodal amplifications and 
will also lead to an enhanced production of photons through the Ua(X) anomaly. 

We have studied the optimal situation in which the order parameter has the largest 
amplitude along the neutral pion direction, with typical amplitudes for the zero mode of 
the neutral pion = 200 — 400 MeV corresponding to energy densities in the initial 
state between 1 — 10 GeV/fm 3 . We have followed the evolution during a "quenched" phase 
transition from an initial state in local thermodynamic equilibrium at a temperature slightly 
above the critical temperature, cooled instantaneously to zero temperature, as well as from 
a non-equilibrium initial condition at Tj = 0. 

We found that the oscillations of the large amplitude neutral pion lead to parametric 
amplification of circularly polarized photons resulting in unstable /c-bands within which the 
photon mode functions grow almost exponentially leading to a distinct distribution of the 
produced photons. The peak of the distribution is correlated with the initial amplitude of 
the neutral pion component of the order parameter and the initial energy density. This is 
a non-perturbative phenomenon as clearly seen in the ratio of the amplitudes of the photon 
distribution functions at the peak between the small amplitude (perturbative) and the large 
amplitude (non-perturbative) regime. 

For initial conditions in which the order parameter has a component only along the n 
direction, we found that the peak in the photon spectrum moves continuously with the 
initial amplitude of the 7r°, from k 300 MeV for (ir°) = 200 MeV corresponding to an 
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initial energy density « 1 GeV/fm 3 , to k ~ 620 MeV for (71*°) = 400 MeV corresponding to 
an initial energy density « 12 GeV/fm 3 during a time scale of 10 fm/c. The pseudoscalar 
nature of the interaction results in a net polarization asymmetry of the final state photons 
which ranges between 0.05 — 0.1 for the above initial conditions. Both the total number of 
photons created during this process and their polarization asymmetry are correlated with a 
non-equilibrium expectation value of the anomaly term FF. 

These novel distributions for photons produced by the "anomalous" decay result in peaks 
in the photon spectrum whose ratio to the normal peak resulting from the "perturbative" 
neutral pion decay at 70 MeV is less than the one percent level. Experimentally this is 
a rather small signal, but perhaps when correlated with the polarization asymmetry, it 
may result in an unambiguous signal for non-equilibrium processes during the chiral phase 
transition. 

We believe that these results are fairly encouraging and justify a deeper study of these 
"anomalous non-equilibrium effects" . The next step in the program is to extend the investi- 
gation of reference [|13[] to understand the effective anomalous vertex out of equilibrium. We 
plan to study this, within the framework of a constituent quark model or alternatively a NJL 
model in which the triangle diagram out of equilibrium should furnish the effective vertex 
for neutral pion decay. After a thorough understanding of the effective tt° — > 27 vertex and 
its translation to an effective mesonic theory, we will include hydrodynamic expansion by 
describing the dynamics in terms of boost invariant variables as in [Q. 
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APPENDIX A: DERIVATION OF NON-EQUILIBRIUM GREEN'S FUNCTIONS: 



In this appendix we present a complete derivation of the non-equilibrium Green's func- 
tions which are needed primarily for the computation of the expectation values of observables, 
such as the photon number density, and FF. 

The Lagrangian density given by eqn. (|3.4j) can be written as 



- A m 
2 



e 2 t(t\ 
S mn n + — ^ J dje mjn 



AVp + total divergence. 



(Al) 



We now define the Fourier transform of the vector potential as, 

d 3 k 



Therefore 
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It also turns out to be convenient to express all our equations in terms of the two operators, 



->mn I zmn 
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k 2 
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(A4) 



where k = y/WW and P mn is the projector onto transverse states. It can be easily verified 
that these operators have the following properties: 
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The action for the electromagnetic sector is then, 
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This is simply a quadratic action for two coupled complex scalar fields. The time dependent 
mixing term, driven by the coherent oscillations of the neutral pion zero mode C{t), ac ts 
like a squeeze parameter. One can decouple the fields, either at the level of the action by 
diagonalizing the quadratic form, or by rewriting the equations of motion in terms of the 
appropriate degrees of freedom. We will adopt the latter method. In any case, the problem 
is completely solved once we obtain all the relevant non-equlibrium two-point functions of 
the theory. The Green's functions Qij for the operator appearing in the quadratic action 
must satisfy: 

V lj g ]k (t,t') = -P lk 5(t-t') (A7) 

where 

°« = jr }■ < A8 > 

The Green's function can be expressed in terms of the two linearly independent operators 
defined in (JA3J) as 

g jk (t, f) = C{t, t')P jk + D{t, t')f jk . (A9) 



Substituting this into eq. (|A7|) , and equating the symmetric and antisymmetric parts of the 
two sides of the equation, we find: 

^ + **-*5^> = -'(«-0 (aio) 

These two differential equations can be decoupled by defining the functions 

G 1 (t,t) = C + iD, G 2 (t,t') = C -iD, (All) 

which then yields 
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G 2 = -5{t-t'). 
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Also, 



Gjk(t, t) 



G\ + G 2 p G\ — G 2 



(A13) 



2 J 2i 

G\ and G 2 now satisfy the simple Green's function equations ( |A12|) and can be readily 
expanded in terms of the mode functions ( |3.8| , |3.9|) as, 



G li2 (t,t')= A>Vu(t)V? a {t') + B>Vl 2 {t)V lt2 {t') Q(t-t') 



(A14) 



+ 



A<V* (t)V 1>2 (t') + B<V h2 (t)V*(t')} Q(t> - t) 



where the mode functions must satisfy the homogeneous equations: 

^ + ev lt(t) - k iLfv^ (ai») 

^ + k*V 2k (t) + *^y^» = 0. (A16) 

The boundary conditions on the Green's functions, which include conditions of continuity 
and the KMS periodicity condition at finite initial temperature, allow all the constants - 
namely, the A's and the B's to be uniquely determined, finally yielding 

M*.0 = 0©(* -0 + ^*(*, t')e(t'-t) (A17) 

with the Wightman functions 

£>(M')= i(A fi ftt)A+ fe (-^,t')) (A18) 
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+ (1 +n k )V 2k {t)V; k {t') + n k V* k (t)V 2k (t')} 



f ~f Ik {(1 + ^)^) W) + n k VUt)V lk (t') 



-(1 + n k )V 2k (t)V* k (t') - n k V* k (t)V 2k (t>)}} 



and 
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Here, and A T represent the transverse components of the gauge field along the forward 
and backward time contours of the closed time path integral. 
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FIGURES 




FIG. 1. n(q,T) = (27r) 3 \q\ d ( N ( T ) ^i ))/^ 1 (j n un its of fm 1 ), vs. q, at time t = 10 fm/c, after 
subtracting the thermal background at T { = 220 MeV, for £(0) = 0.1, £(0) = 0, 0(0) = 0(0) = 0, 
T; = 220 MeV. 
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FIG. 2. (a) £(r) vs. t (in units of fm/c) for f (0) = 1.0, £(0) = 0, 0(0) = 0(0) =0,T t = 220 MeV. 
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FIG. 2. (b) n(q,T) = (2^) 3 |g| <Mll^(gMj ( in units G f fm- 1 ), vs. g, at time t = 10 fm/c for 
the same initial conditions as in Fig. (2. a). 
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2 4 6 8 10 

FIG. 2. (c) h(q,T) = (27r) 3 |g| dN ^J Q (in units of fm -2 ), at the center of the resonant band 
qMp = k = k res = 306 MeV as a function of r (in units of fm/c), for the same initial conditions as 
in Fig.(2.a). 



40 




41 




FIG. 3. (b) n(q,T) = (27r) 3 |g| ^Mll^gMj ( in units G f fm^ 1 ), vs. g, at time t = 10 fm/c for 
the same initial conditions as in Fig. (3. a). 
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FIG. 3. (c) h(q, r) = (2vr) 3 1 e^l — j^'J (in units of fm 2 ), at the center of the resonant band 
qMp = k = k res = 625 MeV as a function of r (in units of fm/c), for the same initial conditions as 
in Fig.(3.a). 
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FIG. 4. (a) E(k = fc res ,r) = (n+ — ra_)/(ra+ + n_) vs. r (in units of fm/c), for fc = 306 MeV 
and for the same initial conditions as in Fig. (2. a). 
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FIG. 4. (b) E(k = k res ,T) = (n+ - n_)/(n+ + vs. r (in units of fm/c), for = 625 MeV 
and for the same initial conditions as in Fig. (3. a). 
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FIG. 5. (b) r(r) = aM F (FF)/8irf 7T (in units of Mj), vs. r (in units of fm/c) for the initial 
conditions of Fig.(3.a). 
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FIG. 6. (a) n(q,r) = (27r) 3 |g| (i n un it s of fm- 1 ), vs. g at time t = 10 fm/c for 
C(0) = 1.0, £(0) = 0, <P(0) = 4>(0) = 0, Tj = 0. 
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FIG. 6. (b) n(q,r) = (2-7r) 3 |g| dN jPJ n (in units of fm~ 2 ), at the center of the resonant band 
qMp = k = k res = 284 MeV as a function of r (in units of fm/c) for the same initial conditions as 
in Fig. (6. a). 
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FIG. 6. (c) E(k = k res ,T) = (n+ - n_)/(n+ + n_) vs. r (in units of fm/c), for k = 284 MeV 
and for the same initial conditions as in Fig. (6. a). 
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